Giant Variations of Cooper-Pair Size in Nanoscale Superconductors 
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The Cooper-pair size (i.e., the BCS coherence length) in low- dimensional superconductors is 
dramatically modified by quantum-size effects. In particular, for nanowires made of conventional 
superconducting materials, we show that the coherence length exhibits size-dependent drops by two- 
three orders of magnitude and reaches values found in high-T c superconductors. This phenomenon 
is surprisingly similar to the well-known BCS-BEC crossover but with one important exception: it 
is driven by the transverse quantization of the electron spectrum rather than by the strength of 
the fermion-fermion interaction. Similar results can be expected for other systems with the same 
structure of the single-particle spectrum, e.g, for superfiuid Fermi gases confined in a quantum-wire 
or quantum-well geometry. 
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Superconductors of ultra-small dimensions possess un- 
usual properties not found in bulk materials. One of them 
is the quantum-size oscillations, first discussed by Blatt 
and Thompson In quasi-lD and -2D superconducting 
systems (nanowires and nanofilms) quantization of the 
transverse electron motion results in single-electron sub- 
bands, i.e., in multiple quantum channels for the super- 
conducting condensate. The proximity of the lower edge 
of a single-electron subband to the Fermi surface leads to 
a size-dependent enhancement of superconducting prop- 
erties, i.e., the superconducting resonance. In particu- 
lar, such resonances are expected to strongly influence 
the critical temperature and critical magnetic field (see, 
e.g., Refs. Q], @, H, 13)- Furthermore, they can lead to a 
remarkable cascade structure of the superconductor-to- 
normal transition Q and can result in the appearance of 
a new type of Andreev states induced by quantum con- 
finement [f|. 

For conventional materials, e.g., Al, Sn or Pb, the su- 
perconducting gap is about 0.1 — 1.0 meV @ and, so, the 
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inter-subband energy spacing ^ (with d the confin- 
ing dimension) becomes of the same order or larger for 
d < 20— 40 nm, where quantum-size oscillations of the su- 
perconducting properties are expected to be significant. 
Several recent experimental results on superconducting 
Pb nanofilms 0, 0] and superconducting aluminum/tin 
nanowires have been attributed to these quantum-size 
effects (see Refs.0, HI and|3, respectively). 

In the present Letter we report an unexpected phe- 
nomenon which is due to quantum-size effects, i.e., gi- 
ant variations of the Cooper-pair size (i.e., the BCS co- 
herence length £o) m low-dimensional superconductors. 
In all previous theoretical studies of superconducting 
nanowires and nanofilms, e.g., modeling phase-slip ef- 
fects in nanowires (see, e.g., Ref. [§), one assumes that 
£o is specified by the same expressions as in bulk. Con- 
trary to this common expectation, our numerical inves- 
tigations of the Bogoliubov-de Gennes equations for 



a clean superconducting quantum wire made of conven- 
tional materials reveal that, depending on the wire width, 
the longitudinal Cooper-pair size varies several orders of 
magnitude, from values of a few micrometers, typical for 
conventional bulk superconductors, to a few nanometers, 
that is usually found in hight-T c materials [10]. This phe- 
nomenon turns out to be very similar to the BCS-BEC 
crossover in superfiuid Fermi gases However, in the 
present case a giant drop in the Cooper-pair size is in- 
duced by quantum-size effects rather than by a change 
in the strength of the fermion-fermion interaction. Our 
results are not only relevant for superconducting quan- 
tum wires but are also applicable to other systems with a 
similar single-particle spectrum, e.g., to ultrathin super- 
conducting metallic nanofilms and ultracold Fermi gases 
confined in a quantum-wire or quantum-well geometry. 

Model and formalism. - We consider a superconduct- 
ing nanocylinder in the clean limit. For our numeri- 
cal calculations we take the material parameters of alu- 
minum, the same as in Refs. 0,0,0]: hw^ = 32.31 meV; 
gN(0) = 0.18 (with AT(0) the bulk DOS); and E F = 

0. 9 eV is the effective Fermi level in the parabolic band 
approximation (for more details, see Ref. 0). Two val- 
ues for the wire diameter are investigated below: d = 
4.22 nm, for which the wire is in the resonance conditions, 

1. e., the bottom of one of the single-electron subbands is 
close to the Fermi level; and d = 4.35 nm, when the wire 
is not influenced by a superconducting resonance. We 
note that our conclusions are not sensible to a particular 
choice of the parameters as long as d < 10 — 15 nm; in 
the wires of larger diameters quantum-size effects play 
less serious role. 

Internal structure of Cooper pairs is described by 
\&(xi,x 2 ) (the Cooper-pair wave function) which is re- 
lated to the anomalous Green's function as 

#(xi,x 2 )=« lim J : (x 1 t 1 ,x 2 t 2 ), (1) 
where Ffati, x 2 t 2 ) = \ {T4>] (xiii)V' i(x2*2)) (for the 
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spin-singlet pairing). In what follows we use cylindrical 
coordinates x = (p, ip, z). The diagonal part of Eq. f2}, 
i.e., at X! = x 2 = x, is directly related to the supercon- 
ducting order parameter A(x) = g^(x, x), where g > 
is the Gor'kov coupling constant. Rotational and trans- 
lational (along the wire) symmetries of the system are 
reflected in \&(xi, x 2 ) which depends on 4> = <p\ — <p% and 
z = z\ — z%. In turn, the order parameter is a function 
of p only. Here, we are interested in the structure of a 
Cooper pair along the wire and, therefore, consider the 
quantity 



*(p, z ) = f, Zl + z \ Pi V, Zl). 



(2) 



The anomalous Green's function can be expressed in 
terms of the eigenstates of the BdG equations (see, e.g., 
Q) which, following the system symmetry, are specified 
by the radial quantum number j, azimuthal quantum 
number to and wavevector k of the quasi-free particle 
motion along the wire. This defines ID subbands labeled 
as (J, to) . For T = the wave function in Eq. ([2]) can be 
written as a sum over such subbands as 



(3) 



with the subband contribution given by 
1 



® jm {p, z) = 



(27T) 



dk u jmk {p)v* mk (p) e lKz . (4) 



Here Uj m k (p) and Vj m k (p) obey the BdG equations writ- 
ten as 
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with Ej m k the quasiparticle energy and 
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where m e is the electron band mass taken equal to the 
free-electron mass. Solutions of Eq. ([5]) are set to be zero 
at the wire boundary (quantum-confinement boundary 
conditions) . 

Resonant subbands, qualitative picture. - The integral 
over k in Eq. (U|) is restricted to the Debye window, i.e., 
\S,jmk\ < hum where £,j m k 1S the subband single-particle 
dispersion £j m k = h 2 k 2 /2m e — pj m with pj m — Ep — 
Ejm the subband chemical potential and Sj m the subband 
lower edge (bottom) energy. Figure [U^a) shows a sketch 
of the single-particle energies for the subbands with j = 
and to = 0, ±1,±2. The dotted horizontal lines in 
Fig. [lja) highlight the Debye window that determines 
the upper k^ m and lower k~ m limits for k in the integral 
in Eq. (gj. 

The bottoms of all single-electron subbands shift in en- 
ergy with changing diameter. A quantum-size supercon- 
ducting resonance occurs when the bottom of a subband 
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FIG. 1: (Color online) (a) Single-electron energies £j m k (mea- 
sured from the Fermi level Ef) versus the wavevector for the 
longitudinal motion k in subbands (j, m) = (0, 0), (0, ±1) and 
(j, ±2). Horizontal dotted lines denote the Debye window, 
(b) Quasiparticle energies Ej m k as function of k for the same 
subbands. 



comes into the Debye window, i.e., when \p.jm\ < fiwn 
and k~ m = 0. In Fig. [IJa) subbands (0,±2) satisfy this 
condition. Below they are referred to as resonant sub- 
bands. Any subband generates a quantum channel for 
the formation of the superconducting condensate. In a 
simplified picture one can utilize Anderson's approximate 
solution of the BdG equations (see, e.g., Ref. [4|), which 
assumes that the spatial dependence of both Uj m k (p) and 
(p) is given by the radial single-electron wave func- 
tion 'djmip) (proportional to the Bessel function of the 
first kind). This leads to E jmk = (£Li. + AD 1/2 , with 
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Ajm the subband energy gap as schematically shown in 
Fig-ETb). 

Within Anderson's approximation, Eq. ((4]) reduces to 
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In most cases the integration limits in Eq. ([6]) can be ex- 
tended to infinity. This yields an exponentially decaying 
function of z, and its characteristic decay length defines 
the subband (channel) BCS coherence length 
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As seen from Eq. (O , {.q" 1 ^ decreases when pj m goes from 
positive to negative values. In the limit pj m 
1 we have the conventional result for the BCS coher- 
ence length, which reads as £q ~ hvj m / Aj mi with 
Vj m — y / 2p,j m /m e the subband Fermi wavevector. At 
resonance, i.e. 



for pj r 



0, we find from Eq. |7]) a 
fi/(m e A jm ) 1/2 . Fi- 



(jm) 



very different expression £g 
nally, when pj m < and \fij m \ ^> A jm , we obtain 
^g- 7 ™ 1 ' ss h/(2m e \pj m \) 1 / 2 1 which decreases with increas- 
ing |/Zj m |. Thus, we have a drop in the BCS coher- 
ence length of the resonant subband(s) and, at the same 
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FIG. 2: (Color online) The resonant wire, d = 4.22 nm: (a) the quasiparticle-energy dispersion Ej m k, the resonant subbands 
are (j,m) = (0, ±7) and (1,±4); (c) the contour plot of \^(p, z)\ (arbitrary units) when accounting only for subbands (0, ±7) 
and (1,±4); (e) the same but for the total Cooper-pair wave function; (g) the longitudinal profile of \^(p, z)\ at p/R = 0.83 
with all relevant quantum channels included (dotted curve) and when only contributions of (0, ±7) and (1, ±4) are taken (solid 
red curve). Panels (b), (d), (e), (f) display the same but for the non-resonant wire at d = 4.35 nm. 



time, such a subband(s) provides a major contribution 
to Eq. §5§ due to a significantly enhanced DOS in the 
Debye window. 

It is surprising that for fij m < 0, Eq. (0 is reduced to 
the expression for the fermion-pair size at the BEC side 
of the BCS-BEC crossover in superfluid Fermi gases (see 
discussion after Eq. (140) in Ref. Illl ). However, in our 
case pij m becomes negative due to quantum-size effects 
instead of a change in the strength of the pair interaction 
(for numerical details, see Fig. [3] and discussion below). 

Numerical solution. - A numerical self-consistent so- 
lution of the BdG equations ([5|) give the results shown 
in Figs. [H and [31 In Fig. a) the quasiparticle ener- 
gies Ej m k are given versus k for d — 4.22 nm. Here 
the bottoms of the two single-electron subbands with 
(j,m) = (0, ±7) and (1,±4) are in the Debye window 
and, so, t>o,±7 and i>i,±4 (recall that Vj m = y / 2pj m /m e ) 
are extremely small. The quasiparticle spectrum for the 
non-resonant wire with d — 4.35 nm is shown in Fig.[2fb). 
The corresponding single-electron spectrum has no reso- 
nant subbands and i>o,±7 and v\ t ±4 are larger by an order 
of magnitude as compared to panel (a). Numerical re- 
sults for the resonant case exhibit a significant increase of 
Ao,±7 and Ai i ±4 and, in turn, lead to enhanced supercon- 
ducting gaps in the other quantum channels. When the 
resonance decays (i.e., due to a change in d), all A jm are 
reduced and approach the bulk value Abuik = 0.25 meV. 

The contour plots in Figs, I^c) and (d), for the res- 



onant and non-resonant wires, respectively, display the 
absolute value of ^(p,z), as calculated from Eq. © 
but with the summation restricted to subbands (0, ±7) 
and (1,±4). Panels (e) [for d = 4.22 nm] and (f) [for 
d = 4.35 nm] show the total quantity |^(p, z)\, where 
we summed over all relevant subbands. For the reso- 
nant case, illustrated by panels (c) and (e), the longitu- 
dinal distribution of electrons is well localized, whereas 
an oscillating and weakly decaying dependence appear 
in panels (d) and (f). This is also clearly seen from 
panels (g) and (h) representing the longitudinal profile 
of \^(p, z)\ at p/R = 0.83 [i.e., the maximum point of 
^{p,z = 0)] for the resonant and non-resonant wires, re- 
spectively. Here the dotted curve gives the total contri- 
bution of all subbands whereas the solid line corresponds 
to a contribution from only (j, to) = (0, ±7) and (1, ±4). 
Thus, the resonant subbands control $>(p,z) for the res- 
onant wire, and the corresponding longitudinal distribu- 
tion of electrons in a Cooper pair is strongly localized. 
At d = 4.35 nm single-electron subbands (j, to) = (0, ±7) 
and (1,±4) are shifted down as compared to their po- 
sitions at d = 4.22 nm. As a result, the resonance dis- 
appears and the relative contribution of the states with 
(j, to) = (0, ±7) and (1, ±4) to the superconducting order 
parameter A(p) = 9(p,z = 0) drops to 40%, see panel 
(h). Nevertheless, the longitudinal decay of the correla- 
tion function z) is still mainly determined by these 
states. 
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FIG. 3: (Color online) BCS-BEC crossover induced by 
quantum-size effects: (a) the subband longitudinal BCS co- 
herence length £g jm ' versus jj,j m as numerically calculated 
from the BdG equations for (j, m) = (0, ±7) (triangles) and 
(1, ±4) (stars), the dashed curve represents Eq. (J7J) with A jm 
replaced by Abuik; (b) the total longitudinal BCS coherence 
length £o versus the nanowire diameter. 



The longitudinal BCS coherence length is defined as 
the decay length of ^(p,z) in the z direction and can 
be calculated through a numerical fit. For both the par- 
tial and total wave functions, such a fitting gives similar 
values: ~ 1 \xyci for the non-resonant wire and s=s 5nm 
for the resonant case. The latter value is almost three 
orders of magnitude less than the BCS coherence length 
in bulk aluminum (« 1.6 /im) and is comparable to the 
one in high-T c superconductors [Io| . Further insight can 
be obtained from Figs. [3£a) and (b). Panel (a) demon- 
strates numerical results for t^™^ as function of fij m for 
(j,m) = (0,±7) (triangles) and (1,±4) (stars). Notice 
that these data are in good agreement with the analyt- 
ical formula of Eq. {?])• When substituting Abuik for 
Aj m in Eq. ([7]), we obtain the dashed curve approach- 

< 0. In this case 



Wel/yml) is n °t sensitive any more to 
Jm . Finally, Fig. [3)Jb) illustrates how the total longitu- 



ing our numerical results for fij 

^' m) « h/(2rn.J,,.^\)V 2 

A 

dinal BCS coherence length £o depends on d. Its value 
is minimal at d — 4.22nm, and the difference between 
the maximum and minimum is roughly two-three orders 
of magnitude. As seen from Fig[3fb), the next super- 
conducting resonance comes into play at d = 4.44 nm. 
Notice that when increasing d, quantum-size oscillations 
of £o are weakened and washed out for d > 30 — 40 nm. 

Conclusions. - We have demonstrated that the lon- 
gitudinal Cooper-pair size in a superconducting metallic 
nanowire undergoes size-dependent giant drops driven by 
the transverse quantization of the electron spectrum. As 



a result, the BCS coherence length of a quantum wire 
made of conventional materials acquires values typical 
for high-T c superconductors. A striking feature of this 
phenomenon is that it is very similar to the BCS-BEC 
crossover found in superfluid Fermi gases. However, there 
is a very important difference: a giant drop of the longi- 
tudinal Cooper-pair size in our quantum superconducting 
wire is a result of the transverse quantization of the elec- 
tron motion while the BCS-BEC crossover in superfluid 
Fermi gases is realized by changing the inter-particle in- 
teraction strength. Notice that the same qualitative be- 
havior can be found for different values of the material 
parameters and a sharp Debye window is not essential. 
Therefore, such a phenomenon is a generic feature that 
will be present in other low-dimensional systems, where 
the condensate is formed via multiple quantum-size chan- 
nels, e.g., in superconducting nanofilms and ultracold 
Fermi gases confined in a quantum- wire or quantum-well 
geometry 



11]. 
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